Velocity selective trapping of atoms in a frequency-modulated standing laser wave 
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The wave function of a moderately cold atom in a stationary near-resonant standing light wave 
delocalizes very fast due to wave packet splitting. However, we show that frequency modulation of 
the field may suppress packet splitting for some atoms having specific velocities in a narrow range. 
These atoms remain localized in a small space for a long time. We propose that in a real experiment 
with cold atomic gas this effect may decrease the velocity distribution of atoms (the field traps the 
atoms with such specific velocities while all other atoms leave the field). 
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INTRODUCTION 

Laser cooling and trapping of atoms and ions is a 
rapidly developing field of modern physics. Cold par- 
ticles in a laser field are a common physical substrate 
used in numerous fundamental and applied issues such as 
Bose-Einstcin condensates, quantum chaos, single-atom 
laser, quantum computer, etc. A significant number of 
methods of atomic cooling in a laser field were devel- 
oped in the recent decades (the Doppler cooling [J Q , 
the Sisyphus cooling 0, 0] , the velocity selective coher- 
ent population trapping (VSCPT) 0, etc. Q). Modern 
sophisticated methods provide temperatures of the order 
of 100 picokclvin 0. 

In this paper we suggest a method of coherent laser 
cooling in the absence of spontaneous emission. Its gen- 
eral ideology is similar to VSCPT in some aspects. The 
field does not cool initially "hot" atoms, it only traps 
the atoms that already have a specific value of veloc- 
ity. However, this velocity is non-zero, and the partic- 
ular trapping mechanism differs from VSCPT radically. 
Our method is not based on " the dark states" . It is 
based on the synchronization between the Landau-Zener 
transitions and the field modulation. When an atom 
moves in a near-resonant standing light wave, two pe- 
riodic optical potentials form in the space [1] . When the 
atom crosses a standing wave node, it may undergo the 
Landau-Zener (LZ) transition between these two poten- 
tials (and the corresponding energy states). Such tran- 
sitions cause splitting of the wave packets Q and rapid 
delocalization of the wave function |10| . In this paper we 
show that frequency modulation of the field may suppress 
the splitting of wave packets for atoms that have veloci- 
ties in the specific narrow range determined by the field 
modulation parameters. We suppose that in a real ex- 
periment, this may significantly decrease the energy dis- 
tribution of moderately cold atoms. This method does 
not pretend to establish any temperature record, how- 
ever, it might be useful in some experiments due to its 



conceptual simplicity. 

In our study the reported effect was initially proposed 
theoretically (semiclassical model) and then confirmed 
numerically (purely quantum model). However, we have 
organized this paper in an alternative order for better un- 
derstanding. First, wc demonstrate the numerical man- 
ifestations of the velocity selective trapping, and then 
explain the effect theoretically. 



EQUATIONS OF MOTION 

Let us consider a two-level atom (with the transition 
frequency Ua and mass rria) moving in a strong standing 
laser wave with the modulated frequency Let us 

assume that the depth of modulation is neglible in com- 
parison with the average value of frequency (a;/[t]) (but 
not with the detuning iOf[t\ — Wq), so we can consider 
the corresponding wave vector fc/ a constant. In absence 
of spontaneous emission (the atomic excited state must 
have long lifetime, or some experimental methods must 
be used to suppress the decoherencc) the atomic motion 
may be described by the Hamiltonian 
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where (T±,2 are the operators of transitions between the 
atomic excited and ground states (the Pauli matrices), 
X and P are the operators of the atomic coordinate and 
momentum, and VL is the Rabi frequency. This Hamil- 
tonian was used in 0-ll|, though for a constant field 
without modulation. 

Let us use the following dimensionless normalized 
quantities: momentum p = P/hkj, time r = Qt, po- 
sition X = kfX, mass m = mai^/hkj: and detuning 
A[r] = (a;/[T] — uja)/^. Let us suppose that the field 
modulation is harmonic. 



A[t] = Ao + AiCos(Cr 



(2) 
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and apply the following conditions: C <C 1, Aq < Ai ^ 
1. Using these approximations we obtain the equations 
for the probability amplitudes to find an atom with the 
normalized momentum p in the excited or ground state, 
a[p, t] and 6[p,t], correspondently: 
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(3) 

Here the dot designates the differentiation with respect 
to T. For every value of p, there is its own pair ([3]). 

Let us choose the values of the parameters and initial 
conditions in order to perform the numerical simulation. 
The average initial atomic momentum (p[0]) will be a 
variable condition for the purpose of this paper. All other 
conditions will be fixed: normalized mass m ~ 10^ (by 
order of magnitude this corresponds to the experiments 
with Cs 12 1 and Rb [l^ atoms, but for a stronger field 
n - 10^-1° Hz), field parameters Aq = -0.02, Ai 
0.047, C = 0.00508, = 0, and the initial form of wave 
packet 
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Therefore, the initial wave packet has a Gaussian form 
with (.^[O]) = and the initial probability to find the 
atom in the excited state 0.5. Here Cp is the standard de- 
viation of the atomic momentum (equal to the half-width 
of the packet by order of magnitude). At r = we fix 
it by the value of ^^[O] = 5\/2. Therefore, in accordance 
with the Heisenberg relation, the standard deviation of 
the initial coordinate is ax[0] = 1/(2<Tj,[0]) = 0.1/\/2 (it 
is much less than the normalized optical wavelength 27r) . 

In numerical experiments, we use these initial condi- 
tions to simulate the system of 8000 equations ([3|) with 
— 1000 < p < lOOO. For larger values of \p\, we put 
a[p, r] = t] ~ due to the energy conservation law. 
Obtaining the solution in the momentum space we per- 
form the Fourier transform and get the wave function in 
the coordinate space in the range of — 47r < x < An (see 
figures in the next section). 



NUMERICAL RESULTS 

In [l3| , the atomic motion was studied in absence of 
field modulation. The following basic modes of motion 
were reported. 

At A = and |A| > 1 the atomic motion is sim- 
ple. Atoms move in constant spatially periodic poten- 
tials. Slow atoms arc trapped in potential wells and fast 
atoms move ballistically through the wave. 



At < |A| ^ 1 the atomic motion is more complex. 
The slowest atoms (|(p[0])| < \j2m) are trapped in po- 
tential wells. Faster atoms {\j2m < \{p[0])\ < 2^/m) 
perform a kind of random walk. Their wave packets split 
each time they cross standing- wave nodes, and this causes 
fast delocalization of the wave functions. The fastest 
atoms (|(p[0])| > 2-ym) move baUistically through the 
wave. Their wave packets split, but all products move in 
the same direction, so the overall delocalization is slow. 
In Fig. [T] we calculate the variance of the atomic po- 
sition cr^ after a relatively long time span of coherent 
evolution t — 5000 as a function of the initial atomic 
momentum (p[0]). For the constant field (solid curve) 
this function shows fast delocalization of all atoms in the 
range of ^/2rn ~ 440 < {p[0]} < 2^/rn ~ 640 (cold atoms 
with velocities of the order of 1 m/s). Local peak at 
(p[0]) ~ 630 is produced by moderately fast atoms hav- 
ing an uncertain scenario of either random walking or 
fiying ballistically. 

Now let us "switch on" the field modulation and see 
the changes. In Fig. [T] the analogous function of cr^ is 
shown with triangles. This function has a more complex 
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FIG. 1: The variance of atomic position cr^ at r = 5000 
as a function of initial atomic momentum (p[0]): curve — 
constant field A(r) = —0.02, triangles — modulated field 
A(r) = -0.02 + 0.047 cos(0.00508r) 



structure. In particular, it has a prominent additional 
minimum at {p[0]) = ptr — 500. These atoms are not 
trapped in potential wells in a strict sense (their energy 
is too high, see the theory in the next sections) , but some 
mechanism significantly suppresses the delocalization of 
their wave functions (note that both functions are shown 
in a logarithmic scale). 

Let us consider the evolution of the corresponding wave 
packets in a coordinate space. In Fig. [2]we show the evo- 
lution of wave functions with (p[0]) = 600 and 500 (other 
parameters are the same as in Fig. [l}. In both cases 
wave packets split. The first splitting occurs near the first 
node, X ~ 1.57 (products overlap at r = 400, but become 
completely independent at r = 800). However, the pro- 
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a: <p[0]> = 600 



b: <p[0]> = 500 
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FIG. 2: Atomic wave packet splitting during quantum evolution (in the coordinate space): (a) fast delocalization of typical 
wave function ((p[0]) = 600), (b) slow delocalization of wave function in the velocity selective trapping mode ((p[0]) — 500). 
Here W[x] is the probability density to find an atom at coordinate x. 



portion of splitting radically differs for {p[0]) — 600 and 
500. In Fig. ^ fission products have similar "weights", 
while in Fig.[5]D they are radically different: a single large 
packet regularly oscillates in the range of — 2 < a; < 2 
"emitting" very small packets in both directions. 



We conclude that the slow delocalization of the wave 
function with {p[0]) ~ 500 is caused by the prominent 
asymmetry of wave packet splitting. Some mechanism 
suppresses the splitting of packets, and the atom is al- 
most completely trapped in the range of — 2 < a; < 2 
(the variance of its position x is even smaller, see Fig.[T]). 
This suppression is significant only for atoms with 490 < 
(-Pp]) ^ 510 (a comparatively narrow momentum and 
velocity range). 



EXPLANATION OF THE EFFECT 

In the previous section we used quantum equations to 
simulate atomic dynamics. In this and the further sec- 
tion, in order to explain the effect of velocity selective 
trapping, let us mention some semiclassical analytical re- 
sults from 0, (obtained for the stationary field). 

In a stationary field with |A| <C 1 the atomic motion 
can be described in terms of two potentials 



A2 



cos^[a;J + — . (5) 



(Fig. [3^1, dashed lines). An atom moves in one of these 
potentials when it is far from the standing wave nodes. 
When an atom crosses the node, the potential may 
change the sign (atom undergoes the Landau-Zencr tun- 
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FIG. 3: (a) Periodic potentials in the space: dashed hne — 
non-resonant potentials solid line — resonant potential 
— cos [a;]; (b) illustration of the trapping condition: the mod- 
ulation of detuning A[r] must be synchronized with atomic 
mechanical motion (A = each time a trapped wave packet 
crosses the standing wave node) 



neling between potentials U^). The probability of tun- 
neling depends on A as exp[— AA^] {A is a combination 
of other parameters). At < |A| <^ 1 the tunneling 
causes splitting of the wave packet (observed in numeri- 
cal experiments). At A = potentials coincide at nodes, 
so the probability of tunneling is equal to 1 and wave 
packets do not split. The correspondent potential takes 
the simplest form U — ±cos[a;] (Fig. solid line). 

What happens, if we "switch on" the field modulation? 
When an atom moves far from the nodes nothing radi- 
cally changes. It moves in a constant potential that does 
not depend much on the value of A. Far from nodes we 
may neglect the term A^/4 in ([5]) and put U ~ ±cos[a;] 
with good accuracy. 

There are two possible scenarios when an atom crosses 
the node (at time r): (1) A[t] 7^ 0, therefore, the packet 
splits significantly; (2) A[t] ~ 0, therefore, the splitting 
is suppressed. 

The first scenario is more typical if the modulation 
is not synchronized with the atomic mechanical motion 
(because most of the time A[t] ^ 0). Second scenario 
may occur sometimes, but does not change the overall 
statistics of the atomic motion. The evolution of the 
wave function shown in Fig. [5^ is typical for moderately 
small detunings |A| ^ 0.01 (both for the stationary and 



the modulated field). 

However, the evolution radically changes if the field 
modulation is synchronized with the atomic mechani- 
cal motion. In particular, it is possible to choose such 
modulation parameters and atomic momentum (the par- 
ticular values are estimated in the next section) that 
A[t] takes zero values each time an atom crosses the 
node. With our parameters such synchronization occurs 
at {p[0]) = ptr ^ 500 (Fig. da). Note that packet split- 
tings are suppressed, but not completely. Slight split- 
tings are caused by the Landau- Zener transitions that 
occur not exactly at a standing wave node, but in its 
small vicinity (when A[t] is small but does not equal to 
zero). 



ESTIMATION OF TRAPPING CONDITIONS 

Let us obtain the analytic relationship between trap- 
ping momentum ptr and field parameters. Trapping oc- 
curs, if A[r] =0 each time atom crosses the nodes of the 
standing wave. In other moments of time |A[r]| ^ 1. 
Therefore, the term A^/4 in ([5]) is always neglible, and 
the trapped atom moves in the effective potential U ~ 
— cos[.t] (we choose the negative sign of [/, because in 
this paper atoms with initial position a;[0] =0 start their 
motion from the potential well). Therefore, the atomic 
center-of-mass motion may be described semiclassically 
by the semiclassical equations of motion 
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i = — , p = — grad[[/] = — sin[x] 
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with the trapping energy 
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being the integral of motion (determined by the initial 
conditions). The trapping energy must be in the range 
of < Etr < 1 (for a:;[0] = 0, this corresponds to \/2m < 
\p[0]\ < 2y/rn). Slower atoms cannot reach the standing 
wave node, and faster atoms move ballistically. 

Let us calculate the atomic traveling time between the 
two successive crossings of nodes in the negative and the 
positive segments of potential — cos x by integrating (|6]) 
(for < Etr < 1) 



r = 2k\/rn, k 
T+ = 2ky/^ ( F 



l+Etr 

TT — I arccos[£'tr] 



(8) 



Here F is the elliptic integral of the first kind. 

In order to synchronize the modulation with the atomic 
mechanical motion the time intervals must be equal to 
time intervals between successive zeros of A[t] (Fig.lJb). 
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Therefore, using ([2]) and ([8]), we get 



c = 
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Ao 
Ai 



cos 



(9) 



These formulae are true for atoms with any initial po- 
sitions (not only a;[0] = used in (|4])). At any value of 
atomic energy in the range of < Etr < 1 (and appropri- 
ate initial momentum) the velocity selective trapping of 
atoms can be achieved with appropriate values of Ao,i, 
C calculated by these formulae. E.g., in order to ob- 
serve trapping at {p[0]) = 500, a;[0] = 0, the field must 
have parameters C = 0.00508, Aq/Ai = -0.4248. We 
use them in numerical experiments, additionally fixing 
Ao = -0.02. 



parameters (see formulae ([8]), ([9])), and only the phase of 
their mechanical oscillations in wells will differ (because 
initial position distribution is random). 

The effect of velocity selective trapping of atoms, being 
theoretically predicted with the semiclassical apparatus, 
has been confirmed by purely quantum numerical mod- 
eling. Therefore, it is not just an artifact of semiclas- 
sical analytics but a real possibility. The drawback of 
this result is that it is obtained in absence of dissipation. 
However, we believe that this is just a quantitative tech- 
nical limitation that may be overcome by an appropriate 
choice of atoms and hi-Q cavities. 

This work has been supported by the Grant of the 
Russian Foundation for Basic Research 12-02-31161. 



CONCLUSIONS 

In this paper we report the effect of velocity selec- 
tive trapping of atoms in a frequency- modulated standing 
laser wave. 

Intensive coherent light produces significant mechan- 
ical action on cold atoms having velocities of the order 
of 1 m/s. There is a wide range of field parameters at 
which atom performs a kind of random walk accompa- 
nied with wave packets splitting and fast delocalization 
of wave function. In this paper we report a specific field 
modulation mode that suppresses wave packet splitting 
for atoms with precisely selected velocities. These atoms 
oscillate in potential wells, and their wave functions arc 
almost completely localized. 

This effect cannot cool atoms in the sense of achieving 
zero velocity, but it can decrease their mechanical energy 
distribution. E.g., if we have a cloud of moderately cold 
atoms having wide position and momentum distribution 
we can switch on the modulated standing wave and wait 
for some time. Most of atoms will leave the wave, and 
only small fraction will be trapped. These trapped atoms 
will have similar mechanical energy determined by field 
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